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A recent study demonstrated the existence of oscillons — extremely long-lived localized config- 
urations that undergo regular oscillations in time — in spontaneously broken SU(2) gauge theory 
with a fundamental Higgs particle whose mass is twice the mass of the gauge bosons. This analysis 
was carried out in a spherically symmetric ansatz invariant under combined spatial and isospin ro- 
tations. We extend this result by considering a numerical simulation of the the full bosonic sector 
of the SU(2) x (7(1) electroweak Standard Model in 3 + 1 dimensions, with no assumption of ro- 
tational symmetry, for a Higgs mass equal to twice the W ± boson mass. Within the limits of this 
numerical simulation, we find that the oscillon solution from the pure SU (2) theory is modified but 
remains stable in the full electroweak theory. The observed oscillon solution contains total energy 
approximately 30 TeV localized in a region of radius approximately 0.05 fm. 

PACS numbers: 11.27. +d ll.15.Ha 12.15.-y 
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Introduction In nonlinear field theories, static soliton 
solutions to the equations of motion have been well stud- 
ied (see for example [H, 0])- However, a much broader 
class of theories contain oscillon solutions, which are lo- 
calized in space but oscillate in time. In some special 
cases, such as the sine- Gordon breather [3( and Q-ball 
[j], conserved charges guarantee the existence of exact, 
periodic solutions. Even in the absence of such guar- 
antees, however, localized solutions have been found in 
many different theories that either live indefinitely or for 
extremely long times compared to the natural timescales 
of the system. 

For scalar theories in one space dimension, oscillon so- 
lutions have been found to remain periodic to all orders 
in a perturbative expansion Q and are never seen to 
decay in numerical simulations [Bj], but can decay after 
extremely long times via nonperturbative effects [|| or 
by coupling to an expanding background 0]. Both </> 4 
theory in two dimensions [a, 0] and the abelian Higgs 
model in one dimension [lOj have also been shown to 
contain oscillon solutions that are not observed to de- 
cay. In </> 4 theory in three dimensions 11, 12, 13|, 14, ll|, 
however, there exist long-lived quasi-periodic solutions 
whose lifetime depends sensitively on the initial condi- 
tions. Similar behavior is present in other scalar theo- 
ries in three dimensions [l 61 ] and in higher dimensions 
[l7j |. Phenomenologically, small_Q-balls were considered 
as dark matter candidates in 
cillons were considered in 122 



and the possible role of 
oscillons in and after inflation was studied in 
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Oscillon-like solutions have also been studied in connec- 
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tion with phase transitions [26j], monopole systems 27j . 
QCD H^i and gravitational systems H3- 

A recent paper [3(| demonstrated numerical evidence 
for an oscillon solution in spontaneously broken SU (2) 



gauge theory with a fundamental Higgs whose mass is 
exactly twice that of the gauge bosons. Current work 
[3~0 ] is investigating an analytic explanation of this mass 
relationship using a small amplitude analysis Hi li 32j |. 
in which the 2 : 1 ratio arises as a resonance condition 
necessary for quadratic nonlinear terms to balance dis- 
persive linear terms in the equations of motion. In this 
analysis, a field of mass m oscillates with amplitude em, 
frequency myl — e 2 , and length scale l/(em). A simi- 
lar mass relation arises in the study of embedded defects 
|33l ]. The field configurations in [3(| were restricted to 



the spherical ansatz 3J], meaning they were invariant 
under combined rotations in space and isospin. Here we 
extend this analysis to a fully three-dimensional spatial 
lattice, eliminating any symmetry assumptions. We in- 
clude the £7(1) hypercharge field, so that we are simu- 
lating the full electroweak sector of the Standard Model 
without fermions. We use the same SU(2) gauge cou- 
pling g and Higgs self-coupling A as in the pure SU(2) 
theory, meaning that the Higgs mass is twice the mass of 
the bosons, and set the £7(1) coupling g' so that the 
mass of the Z° boson matches its observed value. 

While one might expect the oscillon to decay rapidly by 
emitting electromagnetic radiation, it does not. Instead, 
after initially shedding some energy into electromagnetic 
radiation, the system settles into a stable, localized os- 
cillon solution that no longer radiates. Similar behavior 
was observed both when an additional massless scalar 
field was coupled to breathers in one-dimensional <f> the- 
ory and when an additional spherically symmetric mass- 
less scalar field was coupled to oscillons in the spherical 
ansatz, results that provided motivation for this work. 

Continuum Theory We begin from SU (2) x £7(1) elec- 
troweak theory in the continuum, ignoring fermions. We 
follow the standard classical treatment of spontaneously 
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broken nonabelian field theory (see for example [35j|). 
The Lagrangian density is 

(1) 

where the boldface vector notation refers to isovectors. 
Here $ is the Higgs field, a Lorentz scalar carrying U(l) 
hypercharge 1/2 and transforming under the fundamen- 
tal representation of SU(2). The metric signature is 
H . The SU(2) and U(l) field strengths are 



F„ v = d^W v - BJIVp - gW^ x W v , 

Ffj,v = d t _ L B u — d v B^ , 



and the covariant derivatives are given by 



D^Fm = d^F^u - gW^ x F, 



11V j 



J", 
J v , 
2X(v 2 



|$| 2 )$. 



(2) 



(3) 



where r represents the weak isospin Pauli matrices. We 
obtain the equations of motion 



(4) 



Lattice Theory We use the standard Wilsonian ap- 
proach [3^| for implementing gauge fields on the lattice 
(for a review se e 13710 ■ adapted to Minkowski space evo- 
lution as in 39ll40|| . (The details of the discretization 
have been modified slightly for the present application.) 
The U(l) and SU{2) gauge fields live on the links of the 
lattice and the fundamental Higgs field lives at the lattice 
sites. The lattice spacing is Aa;, and we determine the 
values of the fields at time t+ = t + At based on their 
values at times t and t_ = t — At. We associate with 
the link emanating from lattice site p in the positive j th 
direction the Wilson line 



JJP _ e ig'B p ] Ax/2 e igWf-TAx/2 



(9) 



We define the Wilson line for the link emanating from lat- 
tice site p in the negative j th direction to be the adjoint 
of the corresponding Wilson line emanating in the posi- 
tive direction from the neighboring site, U P j = {U p ~ 3 )\ 
where the notation p ± j indicates the adjacent lattice 
site to p, displaced from p in direction ±7. At the edges 
of the lattice we use periodic boundary conditions. We 
define the logarithm of 2 x 2 matrices of this form as 



(10) 



where the gauge currents are 



J v = </Im(A,$)t$, 
J v = g Im(L> l/ $) t r$. 



(5) 



We work in the gauge B = 0, Wo = 0. With this 
choice, the covariant time derivatives become ordinary 
derivatives and we can apply a Hamiltonian formalism. 
The energy density is 



JE 



|2 



k>j 



M 2 + £ (D^iD^ + X^-v 



Fkj ■ Fkj 

2\2 



(6) 



whose integral over space is conserved. Here dot indicates 
time derivative. From the equations for Bo and Wo, we 
obtain the Gauss's Law constraints, 
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E 

j=x,y,z 



D 3 W 3 



Jo 
Jo 



0, 



(7) 



where the charge densities are 



Jo = g'Im&<&, 
Jo = 3 Im $V<]> . 



(8) 



These constraints remain true at all times, at all points 
in space, assuming they are obeyed by the initial value 
data. 



and note that In AB 7^ In A + In B when the matrices do 
not commute. 

The equation of motion for the Higgs field at site p is 



$ p (t+) = 2$p - $ p (i_) + At 2 <t> p . 



(11) 



where 



* P = E + 2A(, 2 -|<F| 2 )<F, 

j=±x,±y,±z 

(12) 

and all fields are evaluated at time t unless otherwise 
indicated. For the gauge fields, we have 



£/?(*+) = exp 



Ax 



lnUfUf(t^ + [^(g'Jf + 9 Jf.T) 



^ hiJ7£,. .,s +ln?7£,. .,,\ 

Q(j,3') o<j,-j') j ^2 



Ax 2 



Uf, (13) 



where U* UJI) = U P U P + J U P + 1+3 U p + J and 



J = g'hn 



J p = g Im 



($p)t[/P$P+j 

Ax ' 
($p)t T [/P$P+j 



(14) 



3 Ax 
are the gauge currents. The energy density at p is then 



E 



1 1 1 exp (lnU p (t + )-lnU p (t_))\ 

2 (2Ai) 2 
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\U P 



-y 

2 ^ Ax 2 

r>j 

AAt 2 



Ax 2 

hAf|$ p | 2 - w 2 ) 



(15) 



whose integral over the whole lattice is conserved. Here 
we have defined 



|Tr lnUjI 



Tr-rln[/, 



P \ t 



(Ttt\uU p ) 



ff ' 2 A:r 2 5 2 Ax 2 
|S?| 2 + Wf -WJ 



(16) 



for any U(2) link matrix. 

At every lattice point, Gauss's Law, 



E 



lnUf(t + )U?(t-)i +lnUZ j (t + )UL j (t-)i 
iAx 2 At 

-(g'jg + gJ$-T)=0, (17) 



is also maintained throughout the evolution, where the 
charge densities are given by 



Jo 
Jo 



g' Im 



g Im 



<f> p (t+) - 3> p (tJ 

2A~t 
$ p (t+) - 
2At 



t 



(18) 



Numerical Simulation The initial conditions for the 
simulation are obtained starting from an approximate 
functional fit to the oscillon solutions found in numerical 
simulations of the SU (2)-Higgs theory in the spherical 
ansatz 30]. This result provides the initial data for the 
W and <E> fields, and the initial B field is chosen to van- 
ish. In order to guarantee that the initial configuration 
continues to obey Gauss's Law in the full SU{2) x U(l) 
theory, we generate the fit at a point in the cycle where 
the time derivatives are smallest, and impose the require- 
ment that all time derivatives vanish for our initial data. 1 

The initial conditions are of the spherical ansatz form 



t ■ Wi 



1 



i \ - - a ( r J t) i - - \ 

ai(r, t)T ■ xxi H (rj — t ■ xxi) 



1 Even for pure SU (2) theory, an approximate fit with nonvanish- 
ing time derivatives will no longer obey Gauss's Law. In that 
case, however, one can restore Gauss's Law by adjusting 
slightly via an SU(2) transformation at each point, 



<■(*+) 



*oId(*+) 



<t> 



with 



exp 



E 



lnUf(t+WW +]nC/*.(t+)(r£,)t 



S 2 Aa: 2 |$ i d (t + )ll$l/2 



(19) 



(20) 



In the St/ (2) xU(l) theory, this procedure is no longer possible 
because <J> carries both charges, so it would be necessary to adjust 
the U(l) field as well, which cannot be done locally. 



(x x r)t 



$ = - [u(r, t) — iv(r, t) r • x] ( !? 
9 V 1 



(21) 



where x = x/r and r is the distance from the ori- 
gin. The field definitions have been chosen so that the 
reduced fields match those used in [3d|, even though 
the conventions for the three-dimensional theory used 
here are slightly different. We work in units where 
v = 1/ v2- Since we are dealing with classical dynam- 
ics, we can rescale the fields to fix the SU(2) coupling 
constant at g = y2, so that the W ± mass is then 
mw = gv/V2 = l/v2- With this rescaling, we must 
also introduce an overall factor of g 1 1 g%/ multiplying the 
total energy, where gw = 0.634 is the true weak coupling 
constant. We choose A = 1, so that the Higgs mass is 
twice the W^ 1 mass, m# = 2v\f\ = \/2. These choices 
of units and normalization agree with |30l] . e xcept the 
quantity called v here is smaller than v in [301 ] by a fac- 
tor of \[2. Finally, we fix g' — 0.773, so that the ratio 
g' jg matches the observed value and the Z° boson has 
the correct mass. 

In these units, we consider initial configurations 



oi(r) 
fi(r) 

v{r) 

a(r) 



e(0.117e + 0.016er) (sech2er) 
1 -0.138esech " 



1/8 



cr 



6.75 ' 



0.017ersech — , 

5_. 

0.117e 2 rsech 
0, 



cr 



(22) 



where the adjustable parameter e allows us to include 
a combined rescaling of the fields' amplitudes and r- 
dependence, as is commonly used in a small amplitude 
analysis 0, [2^, [12] • While e = 1 gives an approximation 
to the spherical ansatz solution of [30], a slightly larger 
value appears to be necessary for the configuration to set- 
tle into a stable solution in the full SU(2) x U(l) model. 
The first term in parentheses in the definition of a\(r) 
is scaled with an additional e so that it matches the co- 
efficient of a, ensuring that a, a\ — a/r, j/r, and v all 
vanish as r — > 0, as required for regularity of the fields at 
the origin. Within the spherical ansatz simulation, these 
initial conditions settle into a long-lived oscillon in the 
pure S'L r (2)-Higgs theory, which we never see decay. The 
U(l) interactions break the grand spin symmetry of the 
spherical ansatz, though the continuum theory still pre- 
serves invariance under combined space and isospin rota- 
tions around the z-axis. The Cartesian lattice provides a 
small breaking of all rotational symmetries. Thus config- 
urations at later times are not constrained to lie within 
this reduced ansatz. (The full three-dimensional simula- 
tion does continue to agree with the spherical ansatz sim- 
ulation for the case when the £7(1) interaction is turned 
off, another check of the numerical calculation.) 
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We start from these initial conditions and let the sys- 
tem evolve for as long as is practical numerically. One 
concern is that the outgoing radiation emitted as the 
configuration settles into the oscillon solution can wrap 
around the periodic boundary conditions, return to the 
region of the oscillon, and potentially destabilize it. How- 
ever, as long as the region in which the oscillon is local- 
ized does not represent a significant fraction of the lattice 
volume, this radiation is sufficiently diffuse that it does 
not affect the oscillon. We use a lattice of size L = 144 
on a side in natural units, which is more than enough 
to satisfy this criterion. For L > fOO, changing the lat- 
tice size simply changes the pattern of noise caused by 
electromagnetic radiation superimposed on the oscillon 
region, but does not affect oscillon properties or stabil- 
ity. We can therefore be certain that there is no coherent 
structure to this unphysical radiation that could possibly 
be necessary for the oscillon's stability; its only possible 
effect is to destabilize the oscillon, which only occurs if 
the radiation is artificially concentrated by a small lattice 
(e.g. of size L < 100). In numerical experiments, these 
destabilization effects are actually much weaker in the 
electroweak model than in pure scalar or SU(2) Higgs- 
gauge models, because in the electroweak model the ra- 
diated energy is almost entirely in the electromagnetic 
field, while the oscillon solution arranges itself to be elec- 
trically neutral. For this reason, it is not necessary to 
use absorptive techniques such as adiabatic damping Q 
(which would have to be adapted to accommodate gauge 
invariance) or an expanding background Q. 

We use lattice spacing Ax = 0.75, though Ax = 0.625 
and Ax = 0.25 were verified to give completely equiva- 
lent results in smaller tests. The time step is At = 0.1. 
Total energy is conserved to a few parts in 10 3 , which is 
appropriate since our algorithm is second-order accurate. 
To check Gauss's Law, we square the left-hand side of Eq. 
(ff"7|) . take the trace, and then take the square root of the 
result. The integral of this quantity over the lattice never 
exceeds 0.025 and shows no upward trend over time, a 
highly nontrivial check on the numerical calculation. It is 
necessary, however, to use double precision to avoid very 
gradual degradation in this result. 2 A run to time 10, 000 
takes roughly 40 hours using 24 parallel processes, each 
running on a 2 GHz Opteron processor core. 

Figure [1] shows the energy in a spherical box of radius 
28 as the fields are evolved from the initial conditions in 
Eq. (j2"2l . When the Higgs mass is twice the W ± mass, 
only a small amount of energy is emitted from the central 
region, with the rest remaining localized for the length of 
the simulation in a stable oscillon solution. If the masses 




- X=0.95 



3,000 12,000 

time 




£ 100 



-X=1 
- A=0.95 



2 Here we have computed the Gauss constraint at time t, which is 
obeyed to order At 2 . The Gauss constraint at time t + At/2 is 
obeyed to machine precision throughout the numerical evolution. 



20,000 30,000 
time 



FIG. 1: Energy in a box of radius 28 as a function of time 
in the natural units of [3C|. One unit of energy is 114 GeV, 
one unit of time is 5.79 x 10 - 27 sec, and one unit of length 
is 1.74 x 10 -18 m. In the top panel, the initial conditions 
are given by Eq. (I22|l with e — 1.15. In the bottom panel, the 
initial conditions are the same except the r z component of W 
in Eq. (|21|l has been set to zero. Two values of the Higgs self- 
coupling A are shown. For A = f , the masses of the Higgs and 
W fields are in the 2 : f ratio needed for oscillon formation and 
the solution remains localized throughout the simulation. A 
transient beat pattern is also visible. For A = 0.95, the mass 
ratio is f.95 : 1. In that case, there is no stable object and 
the solution quickly disperses. 



are not in this ratio, however, the initial configuration 
quickly disperses. The box radius has been chosen to be 
just large enough to enclose essentially all of energy den- 
sity associated with the stable oscillon solution. As a re- 
sult, as the initial conditions settle into the stable oscillon 
solution, we are also able to see a transient "beat" pat- 
tern: the field configurations gradually expand and con- 
tract slightly over many periods, causing a small amount 
of energy to move in and out of the box, accompanied 
by a corresponding modulation of the field amplitudes. 
(When a larger box size is used, the graph of the energy in 
the box flattens out.) Similar beats appear in the SU(2) 
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spherical ansatz oscillon [30|, but in the electroweak os- 
cillon their amplitude decays more rapidly. As in the 
spherical ansatz oscillon, in the electroweak oscillon each 
excited field oscillates at a frequency just below its mass, 
with amplitude of order 0.1 and typical radius of order 
10 in our units. By comparing the total number of cycles 
to the total time, we find Wu = 1.404 for the Higgs field 
components and low — 0.702 for the gauge field compo- 
nents. The primary excitations are in the fields and 
the $ field, with some energy radiated outward in the 
electromagnetic field in a dipole pattern and the Z° field 
largely absent. In contrast, in the spherical ansatz os- 
cillon the W ± and Z° fields must appear symmetrically. 
The electroweak oscillon does remain approximately ax- 
ially symmetric under combined space and isospin rota- 
tions around the z-axis. These results suggest a simple 
modification of the initial configurations in which the r 2 
component of W is set to zero in Eq. (j2Tj) . This modifica- 
tion yields an equivalent final oscillon configuration, with 
the same field amplitudes, localized energy and field fre- 
quencies. However, less energy is shed initially, so there 
is less superimposed noise caused by radiation returning 
from the boundaries. As a result, the beat pattern is also 
more clearly visible. This case is also shown in Fig. [T] 

Conclusions We have seen strong evidence for the ex- 
istence of a long-lived, localized, oscillatory solution to 
the field equations of the bosonic electroweak sector of 
the Standard Model in the case where the Higgs mass is 
twice the M /± mass. Compared to the natural scales of 
the system, this solution has fairly small field amplitudes, 
but because of its large spatial extent it is very mas- 
sive. Such large, coherent objects are well described by 
the classical analysis undertaken here. Quantization of 
the small oscillations around the oscillon solution would 
nonetheless be of interest, perhaps using methods similar 
to those applied to Q-ball oscillons in [41 1. 

Forming oscillons would likely require large energies 
available only in the early universe. In this context, it 
would be very desirable to incorporate fermion couplings, 
which have been ignored here. (Of course, lattice chiral 
fermions introduce significant, but not insurmountable, 
technical complications.) While one might expect the 
oscillon to be destabilized by decay to light fermions, 
in the case of the photon coupling we have seen that 
the analogous decay mechanism is highly suppressed. A 
slow fermion decay mode would be of particular interest 
in baryogencsis, since it would provide a mechanism for 
fermions to be produced out of equilibrium, as is neces- 
sary to avoid washout of particle/ antiparticle asymme- 
try. Or, if the oscillon is extremely long-lived, it could 
provide a dark matter candidate. If such results proved 
compelling, this analysis would suggest a preferred value 
of the Higgs mass. 
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